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We explore how the Peccei-Quinn (PQ) axion parameter space can be constrained by the frequency- 
dependence dimming of radiation from astrophysical objects. To do so we perform accurate calculations of 
photon-axion conversion in the presence of a variable magnetic field. We propose several tests where the PQ 
axion parameter space can be explored with current and future astronomical surveys: the observed spectra of 
isolated neutron stars, occultations of background objects by white dwarfs and neutron stars, the light curves 
of eclipsing binaries containing a white dwarf. We find that the lack of dimming of the light curve of a 
detached eclipsing white dwarf binary recently observed, leads to relevant constraints on the photon-axion 
conversion. Current surveys designed for Earth-like planet searches are well matched to strengthen and improve 
the constraints on the PQ axion using astrophysical objects radiation dimming. 



One of the most attractive solutions to the strong CP 
problem is to introduce the Peccei-Quinn (PQ) symmetry 111. 
This anomalous global symmetry, leads to a pseudo-Nambu- 
Goldstone boson, the axion, whose vacuum expectation value 
permits to solve the problem. Most searches for this pseu- 
doscalar particle are currently being performed in laboratory 
experiments and in astrophysical objects by their generic 
(though model dependent) coupling to two photons ■ 

No PQ axion has been observed yet, which narrows the 
allowed region in parameter space defined by the mass of the 
axion versus its coupling to photons. In the presence of a 
magnetic field, a photon with polarization contained on the 
plane defined by the external magnetic field and the photon 
propagation direction, can be converted into an axion and 
viceversa. This leads, among others, to possible phenomena 
like extra release of energy in stars or appearance of x-rays in 
strong magnetic fields pointing to the Sun. 

It has long been recognized that neutron stars (NS) are 
excellent laboratories to detect axions via the conversion 
effect, because the sensitivity to the axion-photon coupling 
grows with the magnetic field strength dU. Big efforts have 
been focused on axions as dark matter candidates but 
axions may exist and not be the dominant component of the 
dark matter. On the other hand, independently of the fact that 
axion are the dark matter, the PQ axions could leave their 
signature in the photon emission of astrophysical objects such 
as NS and white dwarfs (WD). 

Here, we explore the possibility of detecting PQ axions by 
the frequency-dependent dimming of light of astrophysical 
objects crossing strong magnetic fields. In particular, we 
consider three configurations: ajisolated NS, b) occultation 
of background objects by NS and WD, and c) "eclipsing" 
binary systems containing a WD. We find the last case to be 
the more promising, and also least dependent on astrophysical 
assumptions. Currently on-going surveys driven by the search 
of substellar companions of WD, will provide a large sample 
suitable for this purpose and the first WD-WD eclipsing 
binary has already been observed [jsl). We show that the data 
of [5] akeady place significant bounds on the photon-axion 
conversion. These bounds will be further tested by the ^ 
hundred systems expected to be found by Kepler|l^. 



Photons with polarization in the plane formed by the 
magnetic field and their propagation directions, p||, can be 
converted into axions. The probability P^-^^{L) of a photon 
converting into an axion after travelling a distance L in a 
constant, coherent magnetic field B is given by|0,[ll: 



P-y^4>{L) = sin2(2-0)sin2 
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Here, u! is the photon frequency, M ~ 2TTfaa^^ the coupling 
energy scale {fa being the axion decay constant), a the fine 
structure constant and mofr the effective mass: m^g = — 
ujp I , with wp the plasma frequency of the medium (which 
can be neglected for our purposes) and nia the axion mass. 
Note that AI is related to the axion-photon coupling ga-y-y 
by AI ~ l/ga^y-y. The connection between the PQ mass 
scale and the photon-axion coupling is model-dependent, and 
particle physics constructions ||9| lead to a factor of order 
unity relating the two independent couplings. The function 
q{B) can be a ppro ximated (for 6^1 and b ^ 1) by the 
fitting formula iflOll 



q{B) = 
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where 6 is the magnetic field normalized to the critical QED 
field strength of 4.414 x lO" G. 

However, Eq. ([T]i does not reproduce the correct physics 
in most astrophysical situations, where the magnetic field 
varies along the photon path. For the astrophysical objects of 
our interest, the magnetic field is approximately dipolar (the 
field decreases as distance to the third power) and a correct 
treatment of the inhomogeneous magnetic field is required. 
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Assuming that the magnetic field is constant over distances 
comparable to the photon or axion wavelenghts, the evolution 
of the photon and axion ampUtudes is given by solving the 
Schrodinger-like equation with Hamiltonian given by 
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We obtain the photon-axion conversion by computing numer- 
ically the total evolution operator as the product of evolution 
operators in thin slices (1km width) of constant magnetic field, 



U 



(6) 



The total time, t, from production to detection is subdivided 
in intervals 6tj corresponding to cStj = 1 km. 

We consider three possibilities based on the requirements 
of celestial objects to have large magnetic fields and weak 
dependence on complex astrophysical details: 

a) The peculiar features of the XDINS spectra. These 
isolated NS (e.g., RXJ0720.4-3125, RXJ 1856+5-3754, and 
RXJ 1605. 3+3249 d) have been observed in the eV range 
and the near keV region. A single blackbody spectrum 
fits well the UV region (with temperatures fcr=75,58, 115 
eV respectively), but would yield a too high hard X-ray 
photon flux. With the chosen temperatures, dimming factors 
of roughly 4, 7 and 25 respectively, are needed to fit the 
high-frequency ("X-ray dimmed") data. Simulations of the 
NS magnetic structure predict the photon polarization as a 
function of the frequency, changing from fully perpendicular, 
p±, below few 100 eV to fully p|| at higher frequencies (jsl 
and refs. therein); the frequency range where this transition 
happens is model-dependent. 

b) "Occultation" of background sources by NS and WD. 
If the light of a background source (star or galaxy) passes 
through the magnetic field of a NS or WD, it can be dimmed 
by photon-axion conversion. The relevant radius for a signif- 
icant effect is that of maximum conversion, which is < 10'^ 
(Iff*) km for NS (WD). This is the relevant input for the 
cross section calculation quantifying the expected frequency 
of events. 

c) "Eclipses" in binary systems. This idea was originally 
explored for the double pulsar system detected in 2003 lfl2l . 
Unfortunately the number of binaries with NS is too low 
to be a good candidate to test the conversion into axions. 
We concentrate here in more conspicuous eclipsing binary 
systems where one of the companions is a WD and the other is 
a WD or a main sequence (MS) star The Kepler telescope ||6|, 
aiming at exoplanet transits, should be observing about a 
thousand of eclipsing WD-MS IH], where about 10% of the 
WD have magnetic fields larger than 10^ G. 

Let us begin by discussing the probability of conversion of 
photons into axions produced on the surface of a source. We 
assume a pure beam of p|| photons produced in the presence 
of a large magnetic field; the photon-axion system evolves 
outward towards weaker magnetic filed, where the photons 
are observed. We show in Fig. [T] the dependence of the 
probability conversion into axions for a NS photons and for 
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FIG. 1: Probability conversion of py photons into axions for photons 
exiting a NS dipolar magnetic field as a function of the star's 
magnetic field strength and photon energy. The y-axis denotes the 
NS surface magnetic field and we adopt a typical NS radius of 10km. 
We consider ma=W'^ eV and M = 10^(10'-*) GeV in the top 
(bottom) figure. 



two illustrative M values as a function of surface magnetic 
field and photon frequency. We consider a typical NS radius 
of 10 km; recall that typical surface magnetic fields of NS 
range from 10® G to 10^'^G. In this calculation we account 
for the fact that the photon path is a geodesic (due to the high 
curvature of space time around the NS most light paths are 
not purely radial), but this has negligible effect on the final 
results. In fact most conversion happens at some distance 
for the star's surface (where the ratio of crossing terms to 
the diagonal terms in the Hamiltonian is maximal). Note 
that there are regions with large (90% or higher) suppression 
which are wider for lower coupling energy scales. 

Note the frequency-dependence in Fig. [T] which generally 
leads to a larger dimming effect at high frequencies. White 
regions of the plot correspond to zero conversion efficiency. 
Beside the frequency-dependent conversion, a novel remark- 
able feature is found in our analysis: the fact the maximum 
of the conversion appears at weaker surface magnetic fields 
than the equipartition value. This implies that astrophysical 
systems with weaker magnetic fields than a NS can be suitable 
to explore photon-axion conversion. One excellent candidate 
for harboring such magnetic fields is a WD, about 10% of 
them have magnetic field strengths of about or above lO^G. 
The photon-axion adiabatic conversion is characterized by the 
mixing angle ip in the region of photon production, which is 
small for the weaker magnetic fields. However, for the objects 
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we consider, it is the non-adiabaticity of the photon-axion 
conversion that plays a key role. The physical explanation of 
the net-effect is the non-adiabaticity: the polarization term is 
much smaller than the axion term, the mixing angle ijj changes 
rapidly (as r^'^) along the photon path and the conversion 
happens in the region of maximum non-adiabaticity. 

One important implication of this finding is that, while 
external photons crossing adiabatically a magnetic field do 
not lead to a net conversion into axions, the non-adiabatic 
conversion opens new avenues for the exploration of photon 
dimming. Therefore, we explore the probability of conversion 
of photons into axions crossing the magnetic field of a fore- 
ground star. In this case photons travel through an increasing 
magnetic field at first and later decreasing. The geometry of 
the magnetic field crossed depends on the "impact parameter" 
i.e. the distance from the star's center This will be the case 
of binaries, and/or occultations of background stars by NS or 
WD. We show in Fig. |2] the dependence of the probability 
conversion into axions for two typical values of the coupling 
energy scale M and typical frequency as a function of surface 
magnetic field and radial distance of the transit. This figure 
is qualitatively different from Fig{T] In fact, when a photon is 
generated at the star surface and exits the magnetic field, for a 
given set of parameters ma and M, larger surface magnetic 
field -above a threshold of B - 10^ - 10*G -leads to 
adiabatic photon-axion conversion and therefore to drastic 
reduction of dimming. On the other hand for photons crossing 
the magnetic field of a foreground star, the effect is driven by 
the impact parameter If that is large enough, the photons 
do no cross the non-adiabatic region and no conversion is 
produced, but for moderate surface magnetic fields-S ^ 
lO^G-, outside the star surface there exist a region of non- 
adiabaticity and therefore sizable conversion. 

The two panels of Fig|2] illustrate the sensitivity to the 
axion-photon coupling in this case: the sensitivity to m^ is 
weak for smaller axion masses, while the conversion is lost 
for much larger values. The results shown in Fig. |2]do not 
change with frequency in the optical range. 

In the observation of one of such transits, the axion dim- 
ming effect could be seen as a light-curve dimming. Besides 
statistical errors, the limitation to measure the dimming is the 
accuracy of relative photometry which can easily be better 
than % level. This consideration can be translated into the 
accessible region of the PQ axion parameter space. 

XDINS:We find that a factor of few dimming at high 
frequencies is possible but for values of the magnetic field 
< IQ^G which are severely in disagreement with the values 
obtained from the spin-down analysis (about lO^'^G). Thus 
PQ axions cannot be the explanation for the observed spec- 
trum of XDINS. However, we do find (Fig.l for B - lO^^G) 
that if future observations extend to wavelengths beyond the 
keV range, there could be a noticeable dimming effect even at 
magnetic fields values as the ones measured for XDINS. 

Occultations: Let us now investigate the possibility that 
the effect of photon-axions conversion is observed in the light 
of a background object passing through the influence of NS 
(or WD) magnetic field. We call this effect "occultation". 
A simple estimate of the occultation rate per observed NS 
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FIG. 2: Probability conversion of py photons into axion by crossing 
a star's dipolar magnetic field as a function of surface magnetic field 
strength and distance to the center of the star We adopt a typical WD 
radius of 7000km. The probability regions are weakly dependent on 
the photon energy in the IR-VIS-UV region. We consider 0.3 eV 
photon energy, nia = 10"^ eV and M = 1Q'^(5x 10**) GeV in the 
top (bottom) figure. 



(p') is obtained as follows ll4ll . Consider a NS at distance 
r that passes near the line of sight to a background object 
(which could be a distant galaxy or another galactic star). The 
occultation rate is proportional to the solid angle swept by the 
NS which depends on the projected velocity v and the radius 
from the NS where most of the conversion happens (i?c)- It 
also depends on the number density of background objects on 
the celestial sphere (n^g). Thus we obtain: 
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where in the second line we have rewritten the equation in 
terms of typical numbers for the relevant variables. 

There are iY^g ~ 10^ NS in our galaxy, the total oc- 
cultation rate is therefore obtained integrating over the NS 
distribution. 

We assume lITsll that the NS distribution can be factor- 
ized as function of the distance from the galactic center R 
and height over the galactic plane z, n{R, z)2TrRdRdz — 
N-[^s[nr{R)'^'^RdR][nz{z)dz] . For the R dependence we 
consider two cases. The Hartman exponential model: 
71r{R) = cxp[—R/Rh], where Rh = 5kpc; and the 
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FIG. 3: Lightcurve of an infinitesimal element of the surface of the 
companion when undergoing eclipse behind a WD. We have assumed 
a WD radius of 7000 km, a surface magnetic field of lO^G and the 
two lines are for M = lO'^ (solid) and 3 x lO^GeV (dot-dashed). 
On the X axis we report the projected distance from the center of the 
WD. 



FIG. 4: Region in the PQ rua, Qa-i^ plane excluded by the light 
curve observation of the eclipsing double WD isj] (dark blue). The 
hashed light blue region corresponds to forecasted constraints from 
binary pulsars observations and the black dot denotes the PVLAS 
region. We exclude the possibility of explaining the PVLAS result 
with photon-axion coupling. 



Hartmann model where nii{R) — 25^5- Gxp \^ ^^ 2^"°'' 
with R^, = l.Skpc and i?inax = 3.5kpc; c = 0.204. For 
the z dependence we assume a Gaussian with a = 0.45kpc 
lfl6ll . These distributions must be transformed into spherical 
coordinates (?■, 6, (p) with the Sun at the origin. By integrating 
over the 9 and (p angles we obtain the radial distribution of NS. 
This can then be used to compute the total occultation rate; 



p = / p'n{r)dr = 0.004-0.02 x- 
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1000km 200km/s 100°' yr 
(8) 

depending on the radial NS distribution chosen. 

Therefore the number of expected occultations from NS is 
not significant, even if the whole sky is observed. However, 
we have shown above that maximum conversion takes place 
for surface magnetic fields of order 10^ G. About 10% of 
WD have surface magnetic fields of that order or above. In 
addition, their radius is of the order of 10'' km, which further 
increases the cross section. There are about 5 x 10^" WD in 
the galaxy. So using eq. (15) we find that for WD occultations 
(assuming the same spatial distribution in the galaxy as NS) 
the probability is increased to p = 0.1 — lyr^^, which makes 
the effect much more feasible to be observed. 

Surveys that cover the full sky could find that some of their 
background stars are dimmed when they pass behind a WD 
at a rate of one per few years. No survey currently exists 
that has these characteristics, although in the future Pan-Starrs 
and LSST jT^ should be able to provide few events over 
an operating time period of a decade. It will be interesting, 
nevertheless, to explore if some of such events exist in 
long-running surveys like OGLE. Note that the signature is 
achromatic at optical wavelengths, thus distinguishing it from 
other effects like obscuration. 

WD in "eclipsing" binaries: Another probe to search for 
photon-axion conversion is in binaries where one of the pair is 
a WD or a NS. As it is estimated that about 60% of WD are in 
binaries and WD are much more numerous than NS, we will 



illustrate here the case for WD. 

Any WD in a binary system with orbital plane close to 
the line of sight of the observer would be an excellent can- 
didate for this phenomenon. If the companion's light passes 
within the radius of maximum conversion, Rc, dimming due 
to photon-axion conversion can be observed. Consider an 
infinitesimal element on the surface of the companion star 
and assume it undergoes a full eclipse, although dimming can 
happen as long as the background light passes within Rc- In 
Figure [3] we show the light curve for such an element. We 
assume a typical WD radius of 7000 km, a surface magnetic 
field of lO^G and the two lines are for M = 10^ (solid) and 
3 X lO^GeV (dashed). On the x axis we report the projected 
distance from the center of the WD: this can be mapped into 
time depending on the orbital speed and geometry. In the 
WD- WD binary recently discovered the two WD transit 
in front of each other, and the dimming effect should appear 
in the light curves of both primary and secondary eclipse. 

We fit the phased light curve of the primary eclipse includ- 
ing the photon-axion dimming due to the magnetic field of the 
eclipsing WD. We keep the other light curve parameters fixed 
to the values used in jsll and assume a magnetic field of 10® G 
for the smaller WD. This assumption is motivated by the fact 
that WD in binaries are accreting WD (formed in a common 
envelop) and thus have higher B ifTsIl . In practice, the effect 
illustrated in Fig. 3 is integrated over the projected surface 
of the secondary star. We define a assuming uncorrected 
experimental errors. 

In Fig. |4] we show the resulting 95% exclusion region 
in the ma-ga-y-y parameter space. Using the observed light 
curves, for the first time, we derive bounds on the photon- 
axion conversion by the dimming of light in eclipsing binaries. 
We find that the observations of the WD eclipsing binary, 
subject to the assumption we made on the magnetic field 
strength and the star radius although the region excluded is 
not very sensitive to these parameters, are able to explore and 
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exclude an interesting region of axion parameters: this region 
fully covers the region of direct photon-axion conversion 
explored by laser experiments |fl9| , excludes the possibility of 
explaining the PVLAS results with photon axion conversion 
II20I1 and is competitive with the bounds forecasted in propos- 
als like the observation of high-energy photons in the binary 
pulsar |[l2ll or in high-energy facilities ll2l1l . The bounds 
derived here by analyzing the data of the first WD eclipsing 
binary depend on the radius and the magnetic field of the WD. 



These bounds will be robustly determined by the analysis of 
the data of the large number of eclipsingWD binaries being 
observed by the ongoing Kepler missionja]. 
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